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The hidden ordered state of the frustrated pyrochlore oxide Tb 2 +xTi 2 -cc 07 +^ is possibly one 
of the two electric multipolar, or quadrupolar, states of the effective pseudospin-1/2 Hamiltonian 
derived from crystal-field ground state doublets of non-Kramers Tb^^ ions. These long-range orders 
are antiparallel or parallel alignments of transverse pseudospin components representing electric 
quadrupole moments, which cannot be observed as magnetic Bragg reflections by neutron scattering. 
However pseudospin waves of these states are composite waves of the magnetic-dipole and electric- 
quadrupole moments, and can be partly observed by inelastic magnetic neutron scattering. We 
calculate these spin-quadrupole waves using linear spin-wave theory and discuss previously observed 
low-energy magnetic excitation spectra of a poly crystalline sample with x = 0.005 (Tc = 0.5 K). 


I. INTRODUCTION 

Geometrically frustrated magnets have been actively 
studied in recent yearsi. In particular, pyrochlore 
magnets^ showing spin ice behavior— have interesting 
features such as finite zero-point entropy and emergent 
magnetic monopole excitations^. A quantum spin-liquid 
state is theoretically predicted for certain spin-ice like 
systems^-—, where transverse spin interactions transform 
the classical spin ice into quantum spin liquid. This 
quantum spin ice (QSI), or U(l) quantum spin liquid, is 
characterized by an emergent U(l) gauge field fluctuat¬ 
ing down to T = 0 and by excitations of gapped bosonic 
spinons and gapless photons^iiii^. By changing the in¬ 
teractions of the QSI in some ways the system under¬ 
goes a quantum phase transition to long range ordered 
(LRO) states of transverse spin or pseudospin^, being 
interpreted as Higgs phases^iii. Experimental investi¬ 
gations of the U(l) quantum spin liquid and neighboring 
LRO states have been challenged by several group o^d^d^ . 
However it is difficult to characterize the quantum spin 
liquid states, which preclude standard techniques of ob¬ 
serving magnetic Bragg reflections and magnons. 

Among magnetic pyrochlore oxides^, R 2 Ti 207 (R 
= Dy, Ho) are the well-known classical Ising spin-ice 
examples^. A similar system Tb 2 Ti 207 (TTO) has at¬ 
tracted much attention, because magnetic moments re¬ 
main dynamic with short range correlations down to 
50 mK— . Since TTO has been thought to be close to 
the classical spin ice, the low-temperature dynamical be¬ 
havior of TTO could be attributed to QSV^. Inspired 
by this intriguing idea, many experimental studies of 
TTO have be en p erformed to datei^”— (and references 
in Refs, [ol andl^. However the interpretation of exper¬ 
imental data has been a conundrunt^i^, partly owing to 
strong sample dependenc o^^d^d^ . Among these studies, 
our investigationi^ of polycrystalline Tb 2 +ccTi 2 -cc 07 +y 
showed that a very small change of x induces a quan¬ 
tum phase transition between a spin-liquid state {x < 
—0.0025 = Xc) and a LRO state with a hidden order pa¬ 
rameter {xc < x). It is important to clarify the origin 


of this order parameter, which becomes dynamical in the 
spin-liquid state {x < Xc). 

In this and companio n^^d^ work, we try to reformu¬ 
late the problem of TTO and to reinterpret its puzzling 
experimental data based on the theoretically predicted^^ 
electronic superexchange interactions. A novel ingredient 
of these interactions is the Onoda-type coupling^^ be¬ 
tween neighboring electric quadrupole moments of non- 
Kramers Tb^+ ions. The theory2^ proposes an effec¬ 
tive pseudospin-1/2 Hamiltonian described by the Pauli 
matrices representing both magnetic-dipole and electric- 
quadrupole moments. Depending on the parameters of 
the Hamiltonian there are two electric quadrupole order¬ 
ing phases, which are candidates for the hidden order 
of TTO. These electric quadrupolar orders do not bring 
about observable magnetic Bragg peaks. However, these 
orders can be detected by their elementary excitations 
(inelastic magnetic scattering), and by proper interpre¬ 
tation using a linear spin-wave theory. 

In this paper, starting from the crystal-field (CF) 
ground state doublet of TTO, we account for its single¬ 
site electric quadrupole moments, their LRO, and pseu¬ 
dospin wave excitations in the electric quadrupole LRO. 
A standard linear spin-wave theory predicts that the 
pseudospin wave in the electric quadrupole LRO is, in 
reality, a composite wave of magnetic-dipole and electric- 
quadrupole moments. We discuss this possibility for 
Tb 2 +a:Ti 2 -a^ 07 +y usiug previously observedi^ low-energy 
magnetic excitation spectra of a polycrystalline sample 
with X = 0.005 (Tc = 0.5 K). 

II. CRYSTAL FIELD AND ELECTRIC 
MULTIPOLE MOMENT 

The CF states and inelastic neutron excitation spectra 
of TTO have been investigated by many authors2&“— ; 
readers are referred to Ref. [sO for details. In a low en¬ 
ergy range, there are four CF states: ground doublet 
states and first-excited doublet states at T ^ 16 K. Since 
the interesting temperature range is below 1 K, we ne¬ 
glect the first-excited doublet states and consider only 
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the ground state doublet, for simplicity. 

Among studies of CF, we adopt the CF parameters of 
Ref. [13 (or Ref. [H). The CF ground state doublet of 
TTO can be written by 

I ± 1)d =A\±A)tB\±1)^C\t2)±D\t 5), (1) 

where \m) stands for the \J = 6,m) state within a JLS- 
multiplet^i. The coefiicients2i of Eq. ([1]) are A = 0.9581, 
B = 0.1284, C = 0.1210, D = 0.2256. The local symme¬ 
try axes2^i2S of the crystallographic four sites are 

xo = ^(1,1,2),2 /o = ^(l,l,0),Zo = (2) 

for sites at tn + do with do = |(0, 0, 0), 

= -^{1,1,0),z, = ^{1,1,1) (3) 

for sites at tn + di with di = ^(0,1,1), 

0^2 = ^(1,1, 2), y2 = -^{ 1 , 1,0), ^2 = ;^(1,1,1) (4) 

for sites at tn + c ?2 with d 2 = |(1, 0,1), 

= ^(I, I, 2), ya = ;^(1,1,0), ^3 = 1,1) (5) 

for sites at tn + ds with ds = |(1,1,0), where tn is an 
FCC translation vector. 

In the CF ground state doublet of Eq. o, 
the magnetic-dipole and electric-multipole moment 
operators^ are represented by 2 x 2 matrices: the Pauli 
matrices cr^, cr^, and the unit matrix. The magnetic 
dipole moment operators within | ± 1)d are 

Jx — Jy — 0? 

= (4^2 - 2C^ - bD‘^)(j^ = 3.40c^^ (6) 

which implies that Tb^+ magnetic dipole moments be¬ 
have as Ising-like spins. 

As pointed out in Ref. [HI, for non-Kramers ions in 
the pyrochlore structure including Tb^+ in TTO the CE 
ground doublet states have additionally electric multipole 
moments. These electric multipole moment operators are 
represented by cr^, and the unit matrix. Using the 
explicit form of Eq. m, the electric quadrupole moment 
operators^ within | ± 1)d are expressed by 

i[3j2 - J(J + 1)] = 3A^ - - 15C^ + f 

= 3.05, 

- J2] = (-^B^ + g^lO^c) a* 

= 3.00(7®, 

^[J^Jy + JyJx] = + 9\/T0Ac) ay 

= 3.00(7^, 

^[JzJx + JxJz] = (^-3V^BC - 9^/fAD^ ( 7 ® 

= -8.16(7®, 

^[JyJz + JzJy] = UVmBC + 9 ^Ad\ ay 


Similarly we can show that the electric 16-pole and 
64-pole moment operators^, expressed by the Racah 
operators^ Op^q{J) with p = 4 and 6, respectively (or 
Stevens’s operators), are proportional to or the 

unit matrix within | ± 1 )d. Therefore within the CE 
ground state doublet, pseudospin operators and 
represent the electric multipole moments. A single-site 
CE ground state expressed by 

|tA) = (|l)D,|-l)D)X, (8) 

where y is the pseudospin wave-function, has the largest 
expectation of the magnetic dipole moment |('0|cr^|'0)| = 

1 (and ('0|cr^|'0) = = 0) for y = or y = 

. The other states expressed by 


/cosfe-'^/2\ 
V sin f J 


(9) 


in which 0 is in the range 0 < 0 < tt have finite expecta¬ 
tion values of the electric quadrupole moment operators; 
{'ip\cF^\'ip) 7 ^ 0 and/or {'ip\cry\'ip) ^ 0. These states have 
slightly deformed /-electron charge densities from that 
of the magnetic states with 6> = 0 or 6> = tt. More specifi¬ 
cally, the approximate /-electron charge density^ of the 
state \^) is given by 

(V’Ip(OI^) - (10) 


The angular dependence^ Pe(^) of this equation is 


/9e(f) = n + y; [47r(2p + l)]^/^apYp^g{fydp,g{J), 
p=2,4,6;g 

( 11 ) 

where ( 0 ^ 2 , <^ 4 , tre) = (<^,/^,7) are the Stevens factors^, 
n = 8 is the number of /-electrons, and Yp^q{r) are the 
spherical harmonics. By evaluating {2p\pe{r)\'ip) using 
several spinors of Eq. m, one can show that the defor¬ 
mation of the /-electron charge density is mainly deter¬ 
mined by the electric quadrupole moments. The electric 
16-pole and 64-pole moments have non-negligible contri¬ 
butions to the deformation similarly to the analyses of 
the CE states2&“— . In these meanings, the CE ground 
(psedospin-1/2) states j/i) can be referred to as compos¬ 
ite spin and quadrupole states. 


III. EFFECTIVE PSEUDOSPIN-1/2 
HAMILTONIAN 

The generic form of the effective pseudospin-1/2 
Hamiltonian for non-Kramers CF ground state doublets 
of 4/ magnetic ions in the pyrochlore structure was de¬ 
rived in Ref. [HI by calculating the nearest-neighbor (NN) 
superexchange interaction. This Hamiltonian consists of 
two parts. The first part is the NN magnetic interaction 

-f^m,NN — Ain ^ ^ 

{r,r') 


= 8.16(72^. 


(7) 


( 12 ) 
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which represents the NN classical spin-ice model for 
Jnn > 0- The second part is the NN quadrupolar in¬ 
teraction 

Hq =J„n ^ [2(5((7+cr;:, + cr“cr+,) 

{r,r') 

+ 2g(e2'^--V+a++H.c.)], (13) 

where = (cr^ ± and (a = x^y^z defined 

using the local axes Eqs. (|2])-(|5])) stand for the Pauli ma¬ 
trices of the pseudospin at a site r. The phases (pr,r' 
are (j)ry = 0, —27r/3, and 27r/3 for = (0,3), (1,2), 
{i,i ) = (0,1), (2,3), and = (0,2), (1,3), respec¬ 

tively, where (r, r') = {tn + di, tn' + di>). 

For the magnetic interaction of TTO we probably have 
to include the classical dipolar interaction, i.e., 

Hm = i^m,NN + 

E f Zr • Zr' '^[Zr ' /S.r][Zr' ' Ar] 

\~\K^ \A^ 

{r,r') 

where the summation runs over all pairs of sites, is 
the NN distance, and Ar = r — r' . The parameter D is 
determined by the magnetic moment of the CF ground 
state doublet. We adopt D = 0.29 K, corresponding 
to the experimental value of the magnetic moment 4.6 
As discussed in Refs. |33| and [sj, when the mag¬ 
netic interaction of Eq. (HI represents the dipolar spin 
ice (Jnn > 0), cau be approximated by the NN 

classical spin-ice Hamiltonian^ 




is selected by quantum^ and therma l^^J^ fluctuations for 
PAF, i.e., by an order-by-disorder mechanism. 

In order to calculate elementary excitations in the PAF 
and PF states, we choose one of the pseudospin structures 




and 


(0,(cT*')) (z = 0,3) 
-(0,(<72')) (i = 1,2) (PAF) 

(17) 


«+d.)) = (0, {<y^)) ii = 0,1,2,3) (PF). 

(18) 

We apply the simple linear spin-wave theory, MF-RPA— 
(mean field, random phase approximation), in the same 
way as described in §3.5.2 of Ref. [3l|. In MF-RPA, (cr^) 
of Eqs. (pT|) and ([T8|) is calculated by the MF approx¬ 
imation. For the present purpose, we are interested in 
elementary excitations only at low temperatures, and 
(cr^) = 1 is a good approximation. To obtain dis¬ 
persion relations of pseudospin waves, MF-RPA utilizes 
the generalized susceptibility x(/c,£’) and neutron mag¬ 
netic scattering intensity S{Q^E)^. Useful examples of 
MF-RPA computations including straightforward tech¬ 
nical extensions for pyrochlore structures are described 
Refs. [ 35 I and [36|. General computational treatments of 
MF-RPA are discussed in Refs, and [38|. Following 
these references^i^i^i^, the generalized susceptibility 
is given by 

X{k, E) = [l- x°{E)J{k)]-\°{E). (19) 


alal,, (15) 

{r,r') 

where D^n = \D = 0.48 K. 

In our computations we used an effective pseudospin- 
1/2 Hamiltonian of the form 

= (16) 

We note that this is not very different from the origi¬ 
nal Onoda-t^)e interaction^ (T^nn = 0) and results of 
Refs. [ 25 I and [3 can be approximately used at least in the 
electric quadrupolar phases, in which x^-components of 
the pseudospin (cr^, show LRO and semi-classical the¬ 
oretical treatments are applicable. 


IV. PSEUDOPIN WAVE 

The studies^i^i^ of the effective Hamiltonian of 
Eq. m showed that there are two electric quadrupolar 
states: the PAF state (planar antiferropseudospin) and 
the PF state (planar ferropseudospin) depending on the 
two parameters ((5, q) (see Fig. 7 in Ref. and Fig. 3 
in Ref. 0 for details). In these states, the x^-components 
of the pseudospin show LRO with the modulation vector 
k = 0. It should be noted that this wave vector k = 0 


where k is a wave vector in the first Brillouin zone, x^(F) 
and J{k) denote the single-site generalized-susceptibility 
of the MF Hamiltonian and the Fourier transform of the 
exchange and dipolar coupling constants. The neutron 
magnetic scattering intensity S{Q = G k^E) is given 
by 

S {Q, E) CX ^ “ QpQa) 

p,a 

>< ,(20) 

where only the local ^-component of the pseudospin 
contribute to the scattering. If one assumes that all the 
pseudospin components represent a magnetic dipole mo¬ 
ment vector with an isotropic ^-factor, virtual neutron 
scattering intensity Sy{Q,E) is given by 

S'v (Q, E) OC ^ ~ QpQ<t) 

p,a 

,a' 

where Up]a is the rotation matri x^^i^^ from the local (a) 
frame defined at the sites tn + di to the global (p) frame. 
5'v(Q,F) is useful when displaying dispersion relations 
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FIG. 1. (Color online) Magnetic S(Q,E) (a) and virtual 
Sv{Q,E) (b) of the PAF ordering (Eq. (IT^ ) using interac¬ 
tion parameters Jnn = 1 K, g = 0.85, and (5 = 0. 


FIG. 3. (Color online) Magnetic S{Q,E) (a) and virtual 
Sv{Q,E) (b) of the PF ordering (Eq. ([IH])) using interaction 
parameters Jnn = 1 K, g = 0.8, and S = —0.6. 



EIG. 2. (Color online) Magnetic S{Q,E) (a) and virtual 
Sv{Q,E) (b) of the PAE ordering (Eq. using interac¬ 

tion parameters Jnn = 1 K, g = 0.5, and S = 0.6. 


EIG. 4. (Color online) Magnetic S{Q,E) (a) and virtual 
Sv{Q,E) (b) of the PE ordering (Eq. (HSl) ) using interaction 
parameters Jnn = 1 K, g = 0, and S = —0.6. 


of all pseudospin waves, because the amplitude of the 
electric quadrupole moment are excluded for S{Q,E). 

In Fig. [TJa) we show the inelastic magnetic scatter¬ 
ing intensity S{Q,E) (Eq. ([20])) of the the PAF ordering 
(Eq. (fT7|) ) along several symmetry directions in the FCC 
Brillouin zone using the interaction parameters Jnn = 1 
K, q = 0.85, and (5 = 0 adopted in Ref. S One can see 
two fiat excitation branches in Fig. HJa). We also show 


the virtual Sy{Q,E) (Eq. ([21])) in Eig. Hj^b). This figure 
clearly shows that there are four excitation branches con¬ 
sistent with the k = 0 structure posessing four sites in the 
unit cell. These four pseudospin-wave branches are com¬ 
posite spin (cr^) and quadrupole (cr^) waves. Eigs. (TJa) 
and (b) show that the amplitude of the spin compo¬ 
nents is strong and weak in the two lower-E and the 
two higher-E branches, respectively. Eig. [2] shows the 
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magnetic S{Q,E) and virtual Sy{Q,E) using different 
parameters Jnn = 1 K, q = 0.5, and S = 0.6 in the PAF 
phase (Eq. (pT|) ). The two lower-E excitation branches 
become more dispersive by the finite value of S compared 
to Fig. [H 

In Fig. [3] we show the magnetic S{Q,E) and virtual 
Sy{Q,E) using parameters Jnn = 1 K, g = 0.8, and S = 
—0.6, which are in the PF phase (Eq. ([T8]) ). Compared 
to the PAF cases, the difference between the magnetic 
and virtual S{Q,E) becomes less pronounced. Fig. |4] 
shows the magnetic S{Q,E) and virtual S^^{Q,E) using 
parameters Jnn = 1 K, g = 0, and 6 = —0.6 in the PF 
phase (Eq. (p!8]) b For vanishing g = 0, the two lower-F 
branches are more flattened and merge into almost one 
branch. 


V. MAGNETIC SPECTRA OF 
POLYCRYSTALLINE TB 2 +XTI 2 -XO 7 + 2 / 


Finally, we would like to compare the previously 
observedi^ inelastic magnetic neutron scattering spec¬ 
tra peaked around E = 0.1 meV with the present pseu¬ 
dospin wave calculation. The sample is the polycrys- 
talline Th 2 -\-xT^h-x 07 -^y with x = 0.005 (Tc = 0.5 K)i^. 
The neutron scattering experiment was performed on 
the time-of-flight spectrometer ILL-IN5 operated with 
A = 10 A. Fig. [5fb) shows Q-dependent powder spec¬ 
tra taken at T = 0.1 K. These data should be com¬ 
pared with powder averaging of the magnetic S{Q^E). 
Fig. [5l^a) shows an example of this powder averaged 
S{\Q\,E) choosing the parameters Jnn = 1 K, g = 0.8, 
and (1 = 0, which are in the PAF phase (Eq. ([T7|) b We 
think that these figures show reasonably good agreement 
between the calculation and the observation. In spite 
of using the over-simplified model Hamiltonian for TTO 
and the crude linear-spin-wave theory for the frustrated 
quantum system, essential features of experimental spec¬ 
tra can be reproduced by the approximate calculation. 
The slight Q-dependence and the non-resolution limited 
peak-width, AE ^ (AF)resoiution = 0.01 meV, have been 
one of the puzzling observations of TTO. The present in¬ 
terpretation using the composite spin-quadrupole wave 
can be an answer— 


VI. SUMMARY 


In this study, we try to reformulate the problem 
of Tb 2 +^cTi 2 -a: 07 +y and reinterpret its puzzling ex¬ 


perimental facts based on the theoretically predicted^^ 
pseudospin-1/2 Hamiltonian including the electronic su¬ 
perexchange interaction between electric quadrupole mo¬ 
ments. In this scenario, the hidden order in some TTO 
samples is an electric quadrupolar LRO. Although this 
LRO does not give rise to strong magnetic Bragg scat¬ 
tering, it can be observed by inelastic magnetic neutron 






0.5 n 



0.5 

2 (A-') 



e(A-‘) 


FIG. 5. (Color online) (a) Powder averaged magnetic 
S(\Q\,E) of the PAF ordering (Eq. (ITTF using interaction pa¬ 
rameters Jnn = 1 K, g = 0.8, and (5 = 0. (b) Inelastic neutron 
scattering spectra of polycrystalline Tb 2 +ccTi 2 -cc 07 +y with 
X = 0.005 at T = 0.1 K well below Tc. 


ploy a MF-RPA linear spin-wave theory and compare its 
computation with previously observed low-energy mag¬ 
netic excitation spectra of a polycrystalline sample with 
X = 0.005 (Tc = 0.5 K). Quite intriguingly, the interac¬ 
tion parameters used in Fig. [5fa) are located very close 
to the phase boundary between the PAF and U(l) quan¬ 
tum spin-liquid states^ i^^i^^ . This may possibly imply 
that Th 2 -\-xEi 2 -x 07 -^y samples with x < x^ are in the 
U(l) quantum spin-liquid phase^^. 
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